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Using irreversible thermodynamics we show that current-induced spin transfer torque within a 
magnetic domain implies spin pumping of current within that domain. This has experimental 
implications for samples both with conducting leads and that are electrically isolated. These results 
are obtained by deriving the dynamical equations for two models of non-uniform conducting magnets: 
(1) a generic conducting magnet, with net conduction electron density n and net magnetization M; 
and (2) a two-band magnet, with up and down spins each providing conduction and magnetism. For 
both models, in regions where the equilibrium magnetization is non-uniform, voltage gradients can 
drive adiabatic and non-adiabatic bulk spin torques. Onsager relations then ensure that magnetic 
torques likewise drive adiabatic and non-adiabatic currents - what we call bulk spin pumping. For a 
given amount of adiabatic and non-adiabatic spin torque, the two models yield similar but distinct 
results for the bulk spin pumping, thus distinguishing the two models. As in the recent spin-Berry 
phase study by Barnes and Maekawa, we find that within a domain wall the ratio of the effective 
emf to the magnetic field is approximately given by P(2/is/e), where P is the spin polarization. 
The adiabatic spin torque and spin pumping terms are shown to be dissipative in nature. 

PACS numbers: 72.25.-b, 75.60.Ch, 05.70.Ln 



I. INTRODUCTION 
A. Current-induced Spin Transfer Torque 

Current-induced spin transfer torque at both 
surfaces^! and in bulk is by now a well-established 
phenomenon.^ In surface spin transfer torque, when a 
polarized spin current from a non-magnet crosses the 
interface with a magnet it causes spin motion. In bulk 
spin transfer torque, when a polarized spin current 
crosses a domain wall (where the magnetization varies in 
direction) it causes spin motion. For a uniform magnet, 
there is no current-induced spin torque. 

To be specific, we consider that the magnetization M 
satisfies the equatiorP 

d t M + diQi = -jM xH' + N, (1) 

where Qi is the magnetization flux (or magnetization cur- 
rent density), — jM x H' is the Larmor-like torque, and 
N is the rest of the torque density acting on the magneti- 
zation. We take 7 to be the magnitude of the (negative) 
gyromagnetic ratio. 

H' is the net effective field, which can include an ex- 
ternal field, crystalline anisotropy, the demagnetization 
field, and non-uniform exchange. (We will reserve H 
for a somewhat different quantity, which not only gives 
M x H = in equilibrium, but in fact satisfies H = 
in equilibrium.) H' is given in SI units of T, although 
H' is not the magnetic induction field B, whose units 
also are in T; an H- field in A/m, on multiplication by fiQ, 
becomes an H- field in T. 

We will call 

N' = N- d t Qi (2) 



the non-Larmor-like spin torque. (The prime here is un- 
related to the prime on H'.) When N' contains a term 
that is proportional to the current (or to the gradient 
of the electrochemical potential) one says that there is a 
spin transfer torque. 

For a uniform system, N' will contain only damping 
terms. However, for a nonuniform system N' also con- 
tains terms of the forrrP 

- ZdifiidiM - pfyM x diM], (3) 

where jl the electrochemical potential, ~^dijl has units 
of a velocity, and (3 is dimensionless. The first term is 
called an adiabatic spin transfer torque and the second 
term is called a nonadiabatic spin transfer torque P 

In the literature the current jj usually is written (with 
a suitable conductivity) in place of dijl, but when one 
employs irreversible thermodynamics dijl, a thermody- 
namic driving force that is even under time-reversal, is 
the more natural quantity to employ. Use of dijl un- 
ambiguously leads to the conclusion that, because of 
their time-reversal properties, the adiabatic spin trans- 
fer torque is irreversible (dissipative) and the nonadia- 
batic spin transfer torque is reversible (reactive). This 
also will be seen from a calculation of the volume rate of 
dissipation R, where the equivalent of £ appears, but the 
equivalent of (3 is absent. 

B. Spin Pumping of Current 

Spin pumping of current at surfaces is also a well- 
established phenomenon.^ Here spin dynamics at an in- 
terface transfers a spin-polarized current to an adjacent 
material. The first indication of spin pumping was pro- 
vided by experiments on a thin magnetic film adjacent 
to both a vacuum and an ordinary conductor P The 



present work studies two models for non-uniform con- 
ducting magnets, one a generic conducting magnet that 
is related to but distinct from the s-d modeip and one 
based on a two-band magnet P 

For both models, Onsager relations between transport 
coefficients imply a bulk version of spin pumping, to an 
extent related to the amount of bulk spin transfer torque. 
(This is analogous to how, if temperature gradients can 
cause an electric current, then electrochemical potential 
gradients can cause a heat current, with the size of these 
effects related by Onsager relations.) However, the effects 
are somewhat different for the two models, permitting 
them to be distinguished. 

For the band model, there are two currents (from up 
and down spins) and two effective electrochemical poten- 
tials. Each of these currents can be "spin-pumped" by 
disequilibrium of M. For the generic conducting mag- 
net, there is only one current, but it is taken to be spin- 
polarized. For the generic conducting magnet we study 
the number current density j™ [cf. ( 36 1] . We interpret 
those terms with the form 



Ji 



-L nM1 {diM) ■ H - L' nM diM ■ (M x H) + ... (4) 



as representing the adiabatic and non-adiabatic spin 
pumping terms. Bulk spin-pumping of current requires 
a non-uniform magnetization; it occurs within domains. 

Barnes and Maekawa have recently studied the effec- 
tive electric field associated with the spin-Berry phase 
induced by the domain structure!^ The present results 
are remarkably similar to their, despite the vastly dif- 
ferent methods. In particular, both works find that the 
ratio of the effective emf to the magnetic field is given 
by P(2/is/e), where P is the spin polarization, although 
our result does not appear to be exact. 

The spin pumping terms are related to the spin trans- 
fer torque terms by Onsager relations. (In surface spin 
pumping, there is also a spin transfer torque in propor- 
tion to the spin pumping!^) Thus, not only does the fact 
that spin transfer torque has been observed tell us that 
spin pumping must exist (if we understand the correct 
thermodynamic description of our system) , it also tells us 
how large the spin pumping terms must be. For the two- 
band magnet, analogous terms appear in each of the two 
currents. However, the two models have forms that per- 
mit them to be distinguished. Thus the theory may also 
provide a means to distinguish between the two modelsPJ-' 



C. Summary of Methods 

This paper employs the methods of irreversible ther- 
modynamics. For an introduction to these methods that 
is directed to the magnetism comm unity, s ee RefH"2l We 
summarize the approach as follows.^ 3 * 14 * 15 * 

•The first step is to determine the appropriate ther- 
modynamic variables and the thermodynamic relation for 
the differential de of the energy density as a sum of terms 
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FIG. 1: Experimental geometry to observe spin-pumping of 
current associated with a head-to-head domain wall. An ap- 
plied magnetic field drives the domain wall to the right. A 
spin-pumped current goes to the left. One expects an associ- 
ated voltage pulse when the domain wall crosses either voltage 
lead. 

proportional to the various thermodynamic densities de- 
scribing the system. For an ordinary one-component sys- 
tem these densities are associated with the entropy and 
the particle number. For the generic magnetic conduc- 
tor, these densities are the entropy density s, the density 
of carrier electrons n, and the magnetization M. For a 
two-component non-magnetic system like a semiconduc- 
tor these densities are associated with the entropy and 
the particle numbers (electrons and holes) ™ For the two- 
band magnet, these densities are associated with the en- 
tropy and the particle numbers (up and down electrons, 
with particle densities n^,n^, and the magnetization di- 
rection M. 

•The second step is to require that all of the densities 
satisfy either a conservation law (with unknown fluxes) 
or a source equation (with unknown sources) or both. 

•The third step is to construct the matrix of structure- 
dependent constants (the transport and dissipation coef- 
ficients) relating the sources and fluxes to the thermody- 
namic forces. Essential to the construction of this matrix 
are the time-reversal properties of the fluxes, forces, and 
sources. Finally, the Onsager relations are used to reduce 
the number of independent coefficients. 

The structure of this paper will be, in each section, to 
develop the appropriate parts of each theory in parallel. 



D. Experimental Implications 

The most important prediction of this work, based on 
irreversible thermodynamics, is that bulk spin pumping 
occurs if there is spin transfer torque (and vice- versa); 
moreover, the amount of spin pumping can be determined 
from appropriate Onsager relations. 

For a magnetic wire with head-to-head domains, and 
two voltage leads, if a domain wall is H- field-driven 
past one lead there should be a voltage jump. Figure 
1 shows an experimental geometry corresponding to a 
linear head-to-head domain wall that is driven by an ex- 
ternal field. When the domain wall crosses the voltage 
lead on the right, the voltage difference between the leads 
will measure a pulse that drives current leftward. Like- 
wise, for a magnetic dot with a vortex configuration, and 
two voltage leads, if the vortex structure is TJ-field-driven 
past one lead there should be a voltage jump. 



3 



For Co, data on spin torque indicates that within a do- 
main wall a true E'-field Eq = 1.0 x 10 4 V/m can cause the 
same spin torque as direct application of a magnetic field 
H S t — 4 x 10~ 3 T. From this, in Sect. VII we estimate 
that within a domain wall a true H- field Ho = 0.1 T can 
cause the same spin-pumped current as direct applica- 
tion of an electric field Egp = 350 (V/m), or, across the 
domain wall, can cause the same effect as a voltage dif- 
ference of 7.0 x 10~ 6 V. In terms of the effective emf per 
magnetic field, given by P(2fi B /e), this is 0.7x 10~ 4 V/T. 



II. ENERGY AND ENTROPY 

Energy and entropy are treated the same way in both 
theories. 

Consider the energy density e, which, being conserved, 
has only a flux: 



d t e + dijf = 0, 



(5) 



Here jf is the energy flux density. The intrinsic time- 
reversal signature of e is even, so the intrinsic time- 
reversal signature of jf is odd. 

Now consider the entropy density s. It has both a flux 
jf and a source R/T, where R is the volume rate of heat 
production: 



d t s + dijt 



R 
f 



> 0. 



(G) 



Irreversible thermodynamics considers the time-behavior 
of thermodynamic variables, which have definite signa- 
ture under time-reversal. As a consequence, in an equa- 
tion of motion for that quantity that has both a flux 
and a source, because of the time-derivative, the intrin- 
sic time-reversal signatures of both the flux (here jf) and 
the source (here R/T) are opposite to the intrinsic time- 
reversal signature of then extensive density. Since s is 
even, the intrinsic time-reversal signature of and R/T 
are odd. 



in equilibrium, so that H will serve as a thermodynamic 
force. 

To be specific, consider a uniform ferromagnet in a 
minor loop and in an external field Hq. In the absence 
of anisotropy, in equilibrium it satisfies M = MqHq + 
X-ffo- Let there also be an anisotropy field H an . Further, 
include an a non-uniform exchange term A'(diM) 2 in the 
energy density, which will yield a non-uniform exchange 
field = 2A'V 2 M. We assert that 



H 



5s 
SM 



de 



di 



de 



dM 

= [H + H an + 2A'V 2 M}-[ 
= H' - H 



d(diM) 

M — M M 



X 



lint- (8) 

has the desired form. The first bracket, with three terms, 
constitutes H' , which drives Larmor precession. The last 
bracket, with two terms, represents an internal field H int 
due to exchange. In equilibrium, Hi nt = H' , to yield 
H = in equilibrium. 

Setting H — gives M along H', as desired for no 
Larmor torque. Moreover, it gives M = MqH' + \H ', as 
expected. Even out of equilibrium H int , which is along 
M, does not contribute to M x H, so that M x H' = 
M x H, a substitution we will make when needed. In 
the absence of anisotropy and non- uniform exchange, ^ 
very similar to the form employed by Ref.17, but called 
H*. 

In the s-d model as often used, there is an exchange 
field that couples to the s electrons, giving them a weak 
magnetization. If one wants to incorporate this idea in 
the present framework, then one may consider that M in- 
cludes a contribution from the polarized conduction elec- 
trons. However, because no such specification is made, 
our generic conducting magnet is distinct from the s-d 
model. 



B. Two-band Magnet 



III. THERMODYNAMIC RELATIONS 

A. Generic Conducting Magnet 

We consider a system with magnetization M that is 
basically due to localized electrons, of gyromagnetic ra- 
tio —7, and conduction electrons with electrochemical 
potential jl and number density n. We take the thermo- 
dynamics to be given by 

de = Tds + frdn - H ■ dM. (7) 

Here H includes not only the field H' that causes the 
Larmor torque, but also a uniform exchange field that 
points along M, and thus has no effect on the Larmor 
torque. The exchange field is chose to make H eq = 



The two-band magneiP considers a conducting mag- 
netic system consisting of electrons of charge — e and 
gyromagnetic ratio —7, where 7 = \g\fj,g/H > and 
[iB = eh/2m (for free electrons we take g = —2). We 
assume that the electrons partially occupy two, spin- 
dependent, conduction bands, with number densities rtj 
and ni- For specificity we assume that > nt, so that 
the magnetization will point along the "up" direction, 
determined either by spontaneous symmetry breaking or 
by an external field and anisotropy. The system also has 
an entropy density s. The total number density for the 
conducting electrons is 

n = n ] +n i . (9) 

Moreover, the magnetization is given by 

M = - 7 5, (10) 
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where S is the spin density, of magnitude 

5=(ft/2)(n i -n T ). (11) 

Thus the magnetization has magnitude M = \M\ given 
by 

M = 7 (V2)(«i - «t) = (lffl/V2)(»l - «t). (12) 

With magnetization direction M, we then have 

M = 7 (fy2)(n 1 -n T )M. (13) 

Note that 

dM = 1 {h/2)Md{n l -n ] ) + MdM. (14) 
For the differential of the energy density we take 

de = Tds + ^|dn T + /4«% - M(H ■ dM). (15) 

Here 

^ 4 =MT4±7(V2)(iT-M), (16) 

where 

Mu = Mu- e ^ ( 17 ) 

is the electrochemical potential in terms of the chemical 
potential, with V being the electrical potential. Changes 
in the number density are affected both by the electro- 
chemical potential and H. In equilibrium u* = u* and 

M x H = 0. If we further require that H — in equilib- 
rium, then /ij = //| in equilibrium. 



IV. CONSERVATION LAWS AND EQUATIONS 
OF MOTION 

A. Generic Conducting Magnet 

We take the equations of motion and conservation laws 
for the four variables n and M to be: 

dtn + dij? = 0, (18) 
d t M + d i Q l = - 7 M x H + N. (19) 

Here j™ is the number current density, Qi is the magne- 
tization flux density (t is the real space index), and N (a 
source) is the volume rate of change of magnetization due 
to torques associated with a lack of thermal equilibrium. 
We have employed MxH' = MxH, thus substituting H, 
which in equilibrium is zero, for the precession-causing 
field 



B. Two-band Magnet 

We take the equations of motion and conservation laws 
for the four variables n^, n^, and M to be: 

dtnx+ddv - S, (20) 
dm + dija = S, (21) 
d t M = (7i? + n) x M, (22) 

Here and j^i are the number current densities, S is 
the decay rate for up spins (by charge conservation this 
is compensated by the decay rate —S for down spins), 
and jH and f2 are the Larmor and non-Larmor parts of 
the rotation rate for M. By definition f2 has only two 
components, and is normal to M. (Again we have made 
the allowable substitution, in the torque, of H for H' .) 



C. Two-band Magnet: Implied Equations 

The above equations im ply certain equ ations of mo- 
tion for n of H, M of (fisl and M of H. These equa- 



tions are not necessary, because the previous section is 
self-contained, but they are useful for comparison with 
previous work. 

Continuity equation. Eqs. ([9|, (20), and (21) imply 
that 



d t n+dd? = 0, 



(23) 



With the charge density p — —en and current density 
(charge flux) given by 

P = -en, ji = -ej™ = -e(j Ti + (24) 

the continuity equation is automatically satisfied: 

dtp + dd t = 0. (25) 



Magnitude of Magnetization. Eqs. (12 1, (20), and 



( [21] ) imply that 

d t M + 7(V2)9i(iii - jm) = -25 7 (V2). (26) 
For magnetization along z, this is analogous to dtM z . 



Vector Magnetization. Eqs. (13 1, (20 1, (211, and 



(221 imply that 



d t M+j(h/2)Mdi(jii-j ri ) = --/MxH+nxM-2S-/(h/2)M. 

(27) 

This can be rewritten in the more conventional form 



on setting 



d t M + diQi = -7(M x H) + N (28) 



Qi=j(h/2)M(j li -j n ), (29) 
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where Qi is the magnetization flux density (i is the real 
space index), and 

N = Q x M + 7 (H/2)(j li -j^)d i M-2Sj(H/2)M, (30) 

where N (a source) is the volume rate of change of magne- 
tization due to torques associated with a lack of thermal 
equilibrium. 



A. Generic Conducting Magnet 



From (129]) and ( 30 ) the net non-Larmor spin transfer 



torque of (|2| is given by 

N' = nxM- 7 (h/2)Md i {j li -j ri )-2Sj(h/2)M, (31) 

Thus the only transverse part of N' comes from the f2 x M 
term. 

Once the difference in units are accounted for, Qi above 
is equivalent to Eq. (8) of Ref.O there, both the magneti- 
zation and the magnetization flux densities are measured 
in units of 7, with \g\ = 2, and u is employed for the di- 
rection of th e magnetization. Moreover, the second term 
of N in ( 30 1 is the same as the adiabatic spin torque den- 



sity of Ref. Q30 This adiabatic spin torque is enforced by 
the condition that the magnetization and the spin quan- 
tization axis track with one another. Note that Ref.9 
does not include the adiabatic spin torque density. 



In ( 30 ) , the first and second terms give the transverse 



components in spin space, and the third term gives the 
longitudinal component in spin space. Of course, we have 
yet to determine fl, S, j^, or j^. Below we show that jn 
and ju each have five terms, so that Q a i has ten terms. 
Moreover, has seven terms and S has one term, so that 
N has eighteen terms. 

Note that spin angular momentum is not conserved; 
however, total angular momentum is conserved once one 
accounts for the crystal lattice angular momentum. With 
H' = H + H an + H ex , the angular momentum associ- 
ated with Hq is transferred to the source of Hq (an exter- 
nal magnet or an external current-carrying circuit); the 
angular momentum associated with H an (either lattice 
or dipolar anisotropy) is transferred to the lattice; and 
the angular momentum associated with H ex integrates 
to zero, because it involves the spin system interacting 
with itself. N is associated with the lattice. 



Eqs. (|5|), J6|, fll8| , and placed in the time- 

derivative of (IT]) yields 



< R = 



-3'diT - jfdifi 



HQ, 
cU?4 



N-H. (32) 



Here jf, jf, and Qi are thermodynamic fluxes and diH, 
diT, and dijl are thermodynamic forces; further, N is a 
thermodynamic source, and H is a thermodynamic force. 
In equilibrium all of the thermodynamic forces are zero, 
and thus there is no entropy production. 



Each of the four non-divergence terms in ( 32 ) has a 
clear physical interpretation as a source of heating: the 
first term to thermal conduction, the second to electrical 
conduction, the third to magnetic diffusion (or conduc- 
tion), and the fourth to (local) spin damping. 



B. Two-band Magnet 



Eqs. J6j, and ([20|-(p2| placed in the time- 

derivative of (IT5F yields 

< R = T(d t s + dijf) 

= -9i (it - Tjt - m* it* - vlhi) 

-jtdiT - jading - jadinl 

-(/4 -nl)S+Q- {M x H). (33) 

Here jf, jiji, and are thermodynamic fluxes and diT, 
difj,i, and di/i* are thermodynamic forces; further, S and 

Cl are thermodynamic sources and (n* — fi* ) and M X H 
are thermodynamic forces. 



Each of the five non-divergence terms in ( 33 1 has a 



clear physical interpretation as a source of heating: the 
first term to thermal conduction, the second and third 
to (spin-dependent) electrical conduction, the fourth to 
(local) longitudinal magnetic damping, and the fifth to 
(local) transverse magnetic damping. 



VI. SOURCES AND FLUXES 



A. Generic Conducting Magnet 



V. RATE OF HEAT PRODUCTION 

Irreversible thermodynamics accomplishes its task by 
combining the equations of motion and the thermody- 
namics to obtain an expression for the non-negative 
quantity R as a sum of products of fluxes (or sources) 
and thermodynamic "forces" (or their gradients). 



This model has been treated in the absence of terms 
associated with the conducting electrons.^ 

Energy Flux jf. The energy flux is given by con- 
straining the divergence to be zero (up to an arbitrary 
curl), which leads to 

I I f, ■ /'./," • I' - (I- (34) 

We now must express each flux and source as the sum 
over the suitably weighted "forces" diT, dijl, H , and diH , 
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all of which are zero in equilibrium. The coefficients may 
be constructed from the "order parameters" of the equi- 
librium state, M and d\M. The vector nature of the 
fluxes must be respected (including their properties un- 
der both real space and spin space rotations). 

Entropy Flux jf. The entropy flux, a vector in 
real space whose non-dissipative part is odd under time- 
reversal T, takes the form 

3i = -fdiT - L sn difi- L sQ M ■ diH 

-L sM1 (diM) ■ H - L sM2 (M ■ diM){M ■ H) 
—L' sM diM ■ (M x H). (35) 

There are six terms. The terms with unprimed coef- 
ficients are even under time-reversal T, signifying dis- 
sipation. The term with a primed coefficient is odd 
under T (signifying no dissipation). k is the usual 
thermal conductivity. L sn has units of [diffusion- 
constant-density/ temperature]. L s q, L s mi, and L S M2 
have units of [velocity- length/temperature], or [diffusion- 
constant/temperature]. L' sM has units of [diffusion- 
constant-density / temperature-magnetization] . 

Number Flux j™. The number flux, like the entropy 
flux a vector in real space whose non-dissipative part is 
odd under time-reversal T, takes the form 

3i = -^difi- L ns diT - L nQ M ■ diH 
e* 

-L nM1 {diM) ■ H - L nM2 {M ■ d t M)(M ■ H) 
~L' nM diM ■ (M x H). (36) 

There are six terms. The terms with unprimed coeffi- 
cients are even under time-reversal T, signifying dissipa- 
tion. The terms with a primed coefficient is odd under 
T (signifying no dissipation), a is the usual electrical 
conductivity. The last three terms we will later inter- 
pret as bulk spin pumping terms, and we will relate them 
to the bulk spin torque terms, to be discussed shortly. 
L ns has units of [diffusion-constant-density/ tempera- 
ture]. L n Q, L n Mi, and £„M2 have units of [diffusion- 
constant/energy]. L' nM has units of [diffusion-constant- 
density/ energy-magnetization]. 

Magnetization Flux Q. L . The magnetization flux, 
whose non-dissipative part is even under time-reversal, 
takes the form 

Qi = C\\M{M ■ d l H) - C±M x (M x diH) 
—C'M x diH + L Qs MdiT + L Qn Md t jl 
+L QN1 diM(M ■ H) + L QN2 (M ■ d l M)H 
+L QN3 M(diM ■ H) + L QNA M(M ■ d t M)(M ■ H) 
+L' QN1 diM xH + L' Qm M[(M x d t M) ■ H] 
+L' QN3 (M x diM)(M ■ H) 

+L' QN4 (M ■ diM){M x 3). (37) 

There are thirteen terms. The terms with unprimed co- 
efficients are odd under time-reversal T, signifying dis- 
sipation. The terms with primed coefficients are even 



under T (signifying no dissipation). The terms involv- 
ing H are discussed in Ref. (fl2|l . as is the term involving 
d{I '. Cm and C± have dimensions of a diffusion constant 
divided by a magnetization-squared, and indeed they rep- 
resent longitudinal and transverse diffusion. C has the 
dimensions of [diffusion-constant/magnetization]. Lq„ 
has units of [diffusion-constant/energy]. The term asso- 
ciated with LQ n , when the divergence is taken, will lead 
to a bulk spin torque term. All of the terms in Lqjv and 
L'q N have units of [diffusion-constant/field], with field in 
tesla. Note that, for near-saturation of the magnetiza- 
tion, the terms M ■ diM will be small, so that only Lqni 
and iQjv3 will yield possibly important new terms in the 
damping. 

Spin Torque Density N. The spin torque density, 
whose non-dissipative part is even under time-reversal, 
takes the form 

N = A\\M{M ■ H) - A ± M x (M x H) 

+L Msl d l M8 l T + L Ms2 M(M ■ diM)diT 
+L' Ms (M x diM)diT + L' Mn {M xdiM)diji 

M(M ■ diM)d.Ji 
+L NQ1 M{diM ■ diH) + L NQ2 {M ■ diM)diH 
+L NQ3 diM(M ■ diH) + L NQ4 M(M ■ d t M)(M ■ &H) 
+L' NQ1 diM x d,H + L' NQ2 (M x diM){M ■ &H) 
+L' NQ3 M(M x diM) ■ d t H 

+L' NQ4 (M ■ diM)M x diH. (38) 

There are sixteen terms. The terms with unprimed co- 
efficients are odd under time-reversal T, signifying dis- 
sipation. The terms with primed coefficients are even 
under T (signifying no dissipation). The terms involv- 
ing diH have been discussed in Ref. (fl2)l . as has the 
term involving d{F . The terms in L^nii ^Mnii L' Mn 
we interpret as bulk spin torque terms. and Aj_ 

have units of [l/time-(magnetization) 2 ] for H given in 
A/m; or they must include a factor of 1/(j,q if H 
is given in T. They represent longitudinal and trans- 
verse damping; in terms of the Landau-Lifshitz pa- 
rameter A, one has A± = A/A/. £m s i, and Lms2 
have units of [velocity-length/temperature], or [diffusion- 
constant/temperature]. L' Ms has units of [diffusion- 
constant-density/ temperature-magnetization]. LmuI, 
and Lmu2 have units of [diffusion-constant/energy]. 
L' Mn has units of [diffusion-constant-density/ energy- 
magnetization]. All of the terms in L^q and L' N q have 
units of [diffusion-constant/field], with field in tesla. 

Rate of Entropy Production R. The rate of en- 
tropy production is strictly even under time-reversal. A 
total of forty-one terms can contribute to R when the 
above equations are substituted to find R. The term 
involving C is identically zero. There are six diago- 
nal terms once the longitudinal and transverse parts of 
Qi and N are accounted for. The remaining thirty- four 
terms are cross-terms that occur in pairs. We find that 

R = ^(d t T) 2 + 4(^A) 2 + (L sn + L ns )d z Tdifi 
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+C\\ (M ■ diH) 2 + C±{M x diHf 
+A\\{M ■ H) 2 + Aj_(M x H) 2 
+ (L sQ + Lq s )M ■ diHdiT 
+(L nQ + L Qn )M-d i Hd i fl 
+ {L sM1 + L Msl )d l MjHd l T 

+ {L n Ml + L M nl)diM ■ Hdifl 

+ (L sM2 + L Ms2 )(M ■ diM){M ■ H)d t T 

+ (L nM2 + L Mn2 )(M ■ diM)(M ■ H)diil 

+(L' sM + L' Ms )(Mxd i M)-Hd i T 

+(L' nM + L' Mn )(M x diM) ■ Hdifi 

+ (Lqni + L NQ1 )(M ■ H)(diM ■ d t H) 

+ (L QN2 + L NQ2 )(M ■ diM)(H ■ d z H) 

+ (L Qm + L NQ3 )(M ■ diH){diM ■ H) 

+ {Lqna + L NQ4 )(M ■ O t M)(M ■ H)(M ■ d l H) 

+ (L' QN1 + L / NQ1 )(d l MxH-d l H) 

+ {L' QN2 + L' NQ2 )(M x d t M ■ H)(M ■ ^H) 

+ (L' Qm + L' NQ3 )(M ■ H)(M x d,M ■ diH) 

+(L' Qm + L' NQ4 )(M ■ diM)(M x H ■ diH). (39) 

Recall that H ^ H' , and that H — in equilibrium, 
so that there are no dissipative terms associated with H 
when the magnetization is in equilibrium. 

To ensure that R is invariant under time-reversal the 
non-dissipative cross-terms (which are odd under T) 
must be eliminated. This is done by imposing the condi- 
tions 



L' 



M 



-L' 



A/sj 



NQi 



-Li 



M — 
I 

QNii 



Tl 



Mm 

i = 1,4 



(40) 



In addition, in order to satisfy the Onsager relations 
for dissipative cross-terms (which are even under T) we 
must impose the conditions 

L n s — L ns , L S Q = Lq Si L n Q = LQ n , 
L S M1 = LmsI, L S M2 = Lms2, 
L n Ml — LmuI, L n M2 = ^Mn2! 

Lqnz = L NQi ,i = 1,4. (41) 

We then have 

K 



R 



(diTY 



(d t jiy + 2L sn d z T dip. 



r r '"' ' e 2 
+C\\ (M • diH) 2 + Cj_(M x d.H) 2 

+A\\ (M ■ H) 2 + Aj_(M x H) 2 

+2L sQ M ■ d l Hd l T + 2L nQ M ■ d 1 Hd l fi 

+2L sM xdiM ■ HdiT + 2L nM xdiM ■ Hd^ 

+2L sM 2(M ■ diM)(M ■ H)d t T 

+2L nM2 (M ■ diM)(M ■ H)d z p 

+2L QNl {M ■ H)(diM ■ diH) 

+2L QN2 (M ■ diM)(H ■ diH) 

+2L QN3 (M ■ diH)(diM ■ H) 

+2L QN4 (M ■ d % M){M ■ H)(M ■ d,H) (42) 

Dissipative Nature of Adiabatic Spin Torque. 



are dissipative. We first call attention to the third term, 
due to the thermoelectric effect, involving the product of 
thermodynamic forces c^T dip. Had dip been replaced 
(up to a proportionality constant) by j™, to give the form 
diT j™, which is a product of a thermodynamic force and 
a flux, one would have a term that is apparently odd 
under time-reversal. One therefore might conclude that 
such a term is non-dissipative. Such reasoning would be 
incorrect. In determining the time-reversal properties of 
terms that might contribute to the rate of heating, one 
must work either with products of forces or fluxes; we 
have done the former. 

Similarly, in ( 42 1 the terms in L Mni and Lmu2 , which 



are due both to adiabatic bulk spin torque and adiabatic 
bulk spin pumping terms, produce dissipation. Likewise 
the terms in L n Q, which also produce a (longitudinal) 
bulk spin pumping term and a longitudinal spin flux, 
produce dissipation. 



B. Two-band Magnet 

Energy Flux jf. The energy flux is given by con- 
straining the divergence to be zero (up to an arbitrary 
curl), which leads to 



it = T il - Mtiti - vlhi- 



(43) 



We now must express each flux and source as the sum 
over the suitably weighted "forces" p* — p*, diT, dip*, 

dip*, and M x H, all of which are zero in equilibrium. 
The coefficients may be constructed from the "order pa- 
rameters" of the equilibrium state, M and diM . The 
vector nature of the fluxes must be respected (includ- 
ing their properties under both real space and spin space 
rotations). 

Decay Rate S. The intrinsic signatures under time- 



reversal T of 



and 



are even, so that the intrinsic 



signatures under T of their time-derivatives arc odd. The 
decay rate S is associated with both dtnl and dtn{; hence 
the non-dissipative part of S is odd under T. 
For S, the only possible form is 



S = -a(fi* - fty, 



(44) 



with a a material-dependent constant havings units of 
(J-m 3 -s) _1 . No other form is allowed because S is a 
scalar in both real space and spin space. One might think 
that the "order parameter" M could be multiplied by the 
thermodynamic "force" M x H to obtain a scalar, but 
that dot product is identically zero. Eq. (44 1 is even 



under time-reversal, and therefore is dissipative. 
For small deviations from equilibrium, we have 



-c 



dfi* 
dni 



dul \ 5 n] 

ffn r - 5 -«»il=- 



All the terms in ( 42 ) , being even under time-reversal, 



— = a 



d[i* 
dn^ 



1 

n ' 



dn] 



dn 



(45) 
(46) 
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a result that could have been expected on physical 
grounds. Note that if r-f = r^, then the longitudinal part 
of the magnetization would have a decay rate propor- 
tional to the deviation in the longitudinal magnetization 
itself. However, since under most circumstances local 
electroneutrality enforces Sn^ = —6n^, so SM ~ Sn^, 
this normally will be the case anyway. In this case we 
can think of S as determining T\ processes. 

Entropy Flux j?. The entropy flux, a vector in real 
space whose reversible part is odd under time-reversal T, 
takes the form 



L|| and have the same units. The last two terms in 



31 



T 



d{T - Ls^d^i* - L s idifj,l 



-L sM (M x d,M) ■ (M x H) 
-L' sM diM ■ (M x H). 



(47) 



There are five terms. The terms with unprimed coef- 
ficients are even under time-reversal, signifying dissipa- 
tion. The term with a primed coefficient is odd under 
time-reversal (signifying no dissipation), k is the usual 
thermal conductivity. L s -\ and L s ^, associated with the 
second and third terms, give a well-known electrothermal 
effect, whereby a chemical potential gradient can drive an 
entropy current; they have units of a diffusion constant 
divided by [temperature]. The last two terms, which are 
new, have the same symmetry as corresponding terms in 
the number flux, which we now discuss. They imply that 
spin dynamics can drive an entropy current. We may 
call this spin-pumping of an entropy current. L ^ has 
units of a diffusion constant divided by [magnetization 2 - 
temperature] , and L' ^ has units of a diffusion constant 
divided by [magnetization-temperatu re] . Note that when 
the fifth term is multiplied by d{F in (33 1 for R, unlike R 



the product is odd under T; it thus must be canceled by 
another term (to be discussed below) or it must be zero. 

Number Fluxes and j'u. Each number flux, like 
the entropy flux, is a vector in real space whose reactive, 
or reversible, part is odd under time-reversal T. They 



take the forms [cf. ( 47 1] 



3 H 



a | 

e 



2 di[i\ - L u diT - L tl di^l 



-L ]M (M x diM) ■ (M x H) 
-L\ M diM ■ (M x H), 

-L lM (M x diM) • (A? x H) 
-L\ M diM-(Mx H). 



(48) 



(49) 



Each of these have five terms. The terms with unprimed 
coefficients are even under time-reversal, signifying dissi- 
pation. The term with a primed coefficient is odd under 
time-reversal (signifying no dissipation), a^i gives the 
respective electrical conductivities. L s ^ and L s ±, associ- 
ated with the second and third terms, give a well-known 
electrothermal effect, whereby a chemical potential gra- 
dient can drive an entropy current; they have units of 
a diffusion constant divided by [temperature]; the terms 



both (48 1 and (49 1, which are new, have the same symme- 



try as corresponding terms in the number flux, which we 
now discuss. They imply that spin dynamics can drive 
a current. We call this spin-pumping of current. 
and L.jft have units of a diffusion constant divided by 
[magnctization 2 -energvl , and V - and V - has units of 

a diffusion constant divided by [magnetization-energy]. 
Note that when the fifth terms are respectively multi- 
plied by and c^nj in (33 1 for R, unlike R the product 
is odd under T, and thus these terms must be canceled 
by another cross-term (discussed below) or they must be 
zero. 

Non-equilibrium Rotation Rate VI. The reactive, 
or reversible, part of Vt is odd under T. fl bears the brunt 
of the complexity of the Onsager coefficients. In detail fi 
is given by 



AM x H + L Ms (M x 8iM)diT + L' Ms diMdiT 



-Lmt(M x diM)dni* 
-L ML (M x diM)din\ 



L' M] d i Md l n* 



(50) 



There are seven terms. The terms with unprimed coef- 
ficients are even under T, and thus are dissipative. The 
terms with primed coefficients are odd under T, and thus 
are not dissipative. The first term gives Landau-Lifshitz 
damping. (We believe that any theory based on irre- 
versible thermodynamics that has an energy term vary- 
ing as —HM ■ dM will give a Landau-Lifshitz damping 
term, since then the corresponding thermodynamic force 
is M x H.) The three terms proportional to M x diM are 
even under T, so they are dissipative. They correspond 
to spin torque by electric current and by entropy current. 
Note that there are two types of spin torque, correspond- 
ing to the two types of spins; this is a new result. Both of 
these correspond to what has been called adiabatic spin 
torque. The term driven by entropy current is new. The 
three terms proportional to diM are odd under T, so 
they are non-dissipative. They, too, correspond to spin 
torque by electric current and by entropy current. Again 
there are two types of spin torque, corresponding to the 
two types of spins; this is a new result. Both of these 
correspond to what has been called non-adiabatic spin 
torque. The term driven by entropy current is new. 

Rate of Heat Production R. The rate of entropy 
production is strictly even under time-reversal. We now 



rewrite ( 33 1 in light of ( 43 1 , which eliminates the diver- 



gence term, and in light of the equations for the various 
thermodynamic fluxes and sources. We then have 



R 



A 



A4r + y(W + e ±(%4 

-r — (M x Hf + {diT){di^){L s ^ + L u ) 
+(d i T)(d i 4)(L sl + L ls ) + (difi%)(difj,l)(L n 
+ (M x diM) ■ (M x H) [(L sM + L Ms )diT 

+ (L TM + L M ])diH* + (Aim + L M i)dtH*\] 
+ (diM) ■ (M x H) [{L' sM + L' Ms )diT 



{di»\Y 
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(L'| M + L' M ^)difj^ + (L'^ M + L'^j^di/i*] . (51) Magnetization Source N. This quantity is the sum 



To ensure that R is invariant under time-reversal the 
non-dissipative cross-terms (which are odd under T) 
must be eliminated. This is done by imposing the condi- 
tions 



L' 



-L' 



M i 



L'r 



Ml 



L' 



Ml 



IM- 



(52) 



In addition, in order to satisfy the Onsager relations for 
dissipative cross-terms (which are even under T) we must 
impose the conditions 



L] s — L 



Lms — L s m, 

We then have 



L 



L is — L 



■si- 



L U — L 



M\ 



L 



Ml 



L 



11 ■ 



J.M- 



(53) 



R = - %f + f&T) 2 

A 



+ — (M x Hy + 2L n (d t ^)(d^l) 
+2L sT (9iT)(9 iM p + 2L sl (d i T)(d ifJ r L ) 
+2(M x diM) ■ (M x H) [L sM (d l T) 

+ L TM (9iMt) + L lM{din\)] , 



(54) 



The first five terms are diagonal in the fluxes and source. 
They are always positive. The other six terms, involving 
off-diagonal products of the thermodynamic forces, can 
have either sign, dependent upon the relative directions 
of these gradients. The constraint of positive definitc- 
ness for R then imposes limits the magnitudes of the 
off-diagonal coefficients, which we do not enumerate. 



C. Derived Quantities 



of eighteen terms, which we obtain from (30 1, (48 1, (49 1 



(441, and (50 1. We have 



N= 7 (ft/2) diM (diT(L u -L ls ) 



+%4( e | - L u ) - d lf il(J - L n ) 
-{M x diM) ■ (M x H)(L lM - L TM ) 
-d^ ■ (M x H)(L\ M - L\ M ) 
-( 7 S) M(-l)r(^-/xJ) 

■tfx (\M x H 

+L Ms (M x diM)diT + L' Ms diMdiT 
+L M] {M x diM)di^ + L' m diMdi^ 

+L AU (M x diM)din\ + L' Ml diMdifil 



(57) 



To end this section we note that, by (31), in the ab- 



sence of temperature and chemical potential gradients, 
the non-Larmor spin torque has a transverse part given 
only by the Landau-Lifshitz term. This is in contrast to 
the case of the generic conducting magnet, where two new 
coefficients can contribute to the damping, in a manner 
that depends on the magnetic texture. For the two-band 
magnet there is no such dependence. 



VII. CURRENT-INDUCED SPIN TORQUE AND 
SPIN-PUMPING OF CURRENT 



We now turn to the derived quantities j™, Q a i, and 



Number current density jf. By (23), (48 1, and 



( 49 1 , we have 
3t= -diTiLu+L^-di^Q + L^) 



a *i a l 
d ^l(^ 



-(M x diM) ■ (M x H)(L LM + L ]M ) 
-diM-(MxH)(L\ M + L\ M ). 

Magnetization Flux Q. L . This quantity is the sum 
often terms, which we obtain from (29), (48 1, and (49 1. 
We have 

Q, =j(H/2)m( diT(L^ s — Lis) 

+di^Q - L iT ) - difilQ - L n ) 

— (M x diM) ■ (M x H){L m - L ]M ) 
-d t M ■ (M x H)(L[ M - L\ M )) . (56) 



We are now prepared to discuss both the current- 
induced spin torque and the spin-pumping of current. 
For the two models studied we will restrict ourselves to 
the appropriate terms in the net current density and the 
spin torque. For simplicity we will consider only the sub- 
T^^set of terms that are relevant to spin torque and spin 
pumping, and for that reason we use («) to indicate that 
the appropriate quantities contain certain terms, but are 
(55^iot restricted to those terms. 



Spin Torque and Spin-Pumping Results 

1. Generic Conducting Magnet 



Eq. ( 36 I contains the spin- pumping of number current 
terms 

i? w -L nM1 {diM) ■ H - L nM 2(M ■ 8 t M)(M ■ H) 

-L' nM diM ■ [M x H.) (58) 



The teems involving gradients of the electrochemical po- 
tentials are diffusive in nature. 



The second of these is small if the magnetization is nearly 
saturated, as we will assume. 
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Eq. (38 1 contains the current-driven spin torque terms Eq. (9) of RefJHlhas the form 



N « L' Mn (M x d l M)8 l fi 

+L Mnl d i Md i fi + L Mn2 M(M ■ diMW. (59) 



h at 



HD r - D 



Eq. ( 37 1 contains the current-driven spin torque terms 

Qi w L Qn Mdiji. (60) 

The total non-Larmor spin torque N' thus contains the 
current-driven terms 

N' « L' Mn (M x diM)difl^ 

+(L Mn i - L Qn )d l Md i [i + L M niM(M ■ diM)difl. 

(61) 



For the generic conducting magnet, eqs.(58l and (61) 



provide the basis of our later discussion of the relation- 
ship between spin torque and spin pumping. 



2. Two-band Magnet 

In the two-band magnet, each spin-component has, 
in principle, its own electrochemical potential. In some 
cases such a distinction can be made experimentally. (For 
example, by a suitable combination of electric field and 
of magnetic field gradient it may be possible to produce 
a net spin current but zero net electric current.) For sim- 
plicity, however, let us consider that /i* = /x|. Then, by 
( 55 1 , jf contains the field-driven terms 



j? a -(MxdiM)-(Mx H){L lM + Li M ) 
-diM-(MxH){L' lM + L\ M ). 



(62) 



For the non-Larmor spin torque in the two-band mag- 
net, by (31 1 we need only the transverse ter ms, due to 
Cl x M. If we set //| = /tj = ji, then (50 1 yields the 
current-driven terms 



{L M ] + L Mi ){M x d l M)d l jl 
+(L' m +L' Ml )diMdifl, 



(63) 



so that, by (31 1 
N' 



(L MT + L Ml )M 2 diMdifL 
-{L' m + L' Mi )M 2 {M x diM)difi, (64) 



For the two-band magnet, eqs.(62l and (64 1 provide 



the basis of our later discussion of the relationship be- 
tween spin torque and spin pumping. 



3. Comparison with Previous Work 

Eq. (7) of Refill using a phenomenology based on 
both up and down bands, has the form 



3i 



a 1 + a lEi _D 1± D, din ^ 



(65) 



D ] +D l 



These forms are very similar to what we have derived, in 
that there are three thermodynamic forces in play. For 
RefElthey are the gradient of the voltage, the gradient of 
the density, and the longitudinal gradient of the magne- 
tization. In the two-band magnet they are the gradient 
of the electrochemical potentials, which depend on the 
voltage and on the densities of the up and down spins. 

We now turn to the net torque. The sum of (2) and 
(4) of Rcf.ljhas the form 

N* = -X[M x (M x H)} - v[diM - /3M x ^M], (67) 

where /3 is dimensionless and we expect that (3 <C 1.^ 

Microscopic theories for the adiabatic spin torque give, 
with P the polarization of the current (e.g., 0.6 for 60% 
in the up band)^ 



eM 



(model) . 



(68) 



To make a proper comparison with this form, we must 
replace the current density ji (a flux) by the form it 
takes when driven by the "force" dijl. We also employ a 
less model-dependent form by introducing the constant 
£, with units of a diffusion constant divided by energy, 
and let ji — > (a/e)dijl. Thus we write 



(69) 



so (67 1 becomes 



N* 



-X[Mx(MxH)]+(d i fi[d i M-f3Mxd i M]. (70) 



In this form, which is appropriate to irreversible ther- 
modynamics, the adiabatic spin torque (proportional to 
diM) is odd under time-reversal, opposite the even sig- 
nature of a non-dissipative spin torque. Therefore the 
adiabatic spin torque is dissipative, as can be seen from 
its contribution to the rate of heating R above. On the 
other hand, the non-adiabatic spin torque (proportional 
to M x diM) is even under time-reversal, signifying that 
it is non-dissipative, as can be seen by its absence from 
R. 

For the above model we then have 



_ PauB 
4 ~ e 2 M 



(model). 



(71) 



To close this subsection we note that Ref . |7] uses very 
different methods to show that, at surfaces, spin trans- 
fer torque and spin pumping are related. We also note 
that Ref. [10] uses a Berry-phase and a spin-Berry-phasc 
to predict that, for ferromagnetic conductors, there is 
an effective spin-dependent emf that drives an ordinary 
electric current and a spin emf that drives a spin current. 
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4- Comparison with the two models 



Now define 



Generic Conducting Magnet. Comparison of (70) 



and the generic conducting magnet result (61 1 shows that 



Rst = 



Hst_ 
~E~ 



(78) 



the two versions of the spin torque have the same form if 
(L Mn i - L Qn ) = £, L' Mn = — /3£. (72) 



which is in T-m/V. Then, by (|75j) and (|76j) , eq. (|78| yields 

Nst T e 



Rst 



jME 



' 7 cT 



(79) 



Two-band magnet. Comparison of (70 1 and the two- 



the spin torque have the same form if 

(L MT + L Mi )M 2 = £, {L' m + L' Mi )M 2 = ft. (73) 

Note that the generic ("g") conducting magnet and 
the two-band ("2b") magnet have transport coefficients 
with different dimensionality. By ( 72 I and ( 73 1 , for later 
purposes we write 



band magnet result (64) shows that the two versions of S p m torque terms of (|70| as 



From the "velocity" of ( 68 ) and ( 75 1 , we can write the 



L g » L 2b M 2 » e- 



Estimates 



(74) 



For purposes of estimation we will employ Table 1 of 
Ref. M which for Co gives P = 0.6, M = 1.45 x 10 6 A/m 
and implies that a = erj+crx = 3.4 x 10 7 /fkm. With stan- 
dard values of e = 1.6 x 10~ 19 C and fi B = 9.3 x 10~ 24 A- 
m 2 , we estimate that, for Co, £ » 5.1 X 10 15 m 2 /J-s. (As 
indicated above, this has units of a diffusion constant over 
an energy. Indeed, use of a = ne 2 T/m and M — Pn[iB 
gives £ = r/m. For a bare electron mass m this corre- 
sponds to the somewhat short time of r w 4.5 x 10~ 15 s.) 
Also note that 7 = 1.9 x 10 n /T-s. 

Let (5 be a characteristic domain wall dimension, which 
for purpose of estimation (in Co) we will take to be 
10 nm. Within the domain walls we may take \diM\ w 
2M/8 (2 because the magnetization reverses). In order 
to make estimates, we will neglect the vector nature of 
various quantities, and not consider signs. We will con- 
sider only the adiabatic spin torque, which we believe 
dominates experimentally! 1 ^! 

On neglecting spin-pumping terms a field E yields a 
current density j proportional to the conductivity a via 



j « oE. 



(75) 



For non-uniform magnets there is also a current- 
induced spin torque N ST that is proportional to j. We 
define the equivalent spin torque field Hst via 



H 



ST 



N ST 
7M ' 



(76) 



Fro m (|70[ ), with L g the appropriate Onsager coefficient 
(by (1721), for estimation purposes L g is like £, with units 
of a diffusion constant divided by an energy), we have 



T eM E 



(77) 



Nst — PoE^b 



1 M _ Pam 
eM 6 ~ eS 



E. 



Comparison of (77 1 and (|80| gives the result 
L„ 



Pans 



e 2 M 



(80) 



(81) 



This permits the estimate that, for Co, L g = 5.1 x 
10 15 (m 2 /J-s). Then, again for Co, application of (79 1 
gives Rst = 4 x 10~ 7 (T-m/V). Therefore within a do- 
main wall a true -E-field Eq = 1.0 x 10 4 V/m can cause 
the same torque as direct application of a magnetic field 
Hst = 4 x 10~ 3 T. This merely restates what is already 
known. 



From ( 58 1 , an applied field Hq that produces a torque 

N = jMH (82) 

also produces a spin-pumping driven current within a 
domain wall. With \diM\ « 2M/S, this current is given 

by 



2 M 

] S p = eL g — H Q . (83) 
We define the equivalent spin-pumping field Esp via 



E. 



SP 



With the definition 



Rsp 



JSP 

a 

Esp 
~Hq~ 



we find that 



Rsp = 



JSP 



2eM 



(84) 



(85) 



(86) 



Note that Rsp and Rst have units that are inverse to 
one another. 

For Co this gives R S p = 7.0 x 10 3 (V/T-m). Thus, 
within a domain wall a true -ff-field 7?o=0.1T can cause 
the same current as direct application of an electric field 
Esp — 700(V/m). This is a new prediction. The corre- 
sponding voltage difference across the domain wall is on 
the order of E SP S, or 7.0 x 10~ 6 V. 

Consider the situation depicted in Figure 1, where 
H points rightward and d x M points leftward. With 
L n Mi ~ £ > 0, by (58 1 the number current density points 
rightward, so that (within the domain wall) the electric 
current density points leftward. 
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C. Comparison with Barnes and Maekawa 

For comparison with Ref. 1101 we define a spin- pumping 
emf 



&sp — EspS. 



(87) 



Then the emf divided by the field, on using (851, ( 86 1 
and (81 1, is given by 



S S p = E SP § = 2P VB 
Ho H e 



Ref. [T0]finds a quantity £ s to be given by 2^BH /e, and 
that P£ s drives a current density. On division by H, this 
is precisely (88). We believe that this exact agreement, 



but not the parameter dependence, is accidental. 



VIII. SUMMARY AND DISCUSSION 

Using irreversible thermodynamics we have shown, for 
both a two-band magnet and a generic conducting mag- 
net, that current-induced spin transfer torque within a 
magnetic domain leads, by Onsager relations, to spin 
pumping of current within that domain. For a given 
amount of adiabatic and non-adiabatic spin torque, the 
two models yield similar but distinct results for the bulk 
spin pumping, thus distinguishing the two models. 

This has experimental implications both for samples 
with conducting leads and that are electrically isolated. 
For Co we estimate that within a domain wall a true 
ii-field Ho = 0.1 T can cause the same spin-pumped 



current as direct application of an electric field E$p — 
350 (V/m), or, across the domain wall can cause the same 
effect as a voltage difference of 7.0 x 10~ 6 V. Correspond- 
ingly, the ratio of the effective emf to the field Ho is, for 
Co, about 0.7 x 10~ 4 V/T. 

The similarity between our results and those of Barnes 
and Maekawa is likely not an accident. In the present case 
we have shown that the "off-diagonal" current-induced 
( "adiabatic" ) spin transfer torque implies a similar "off- 
diagonal" ("adiabatic") spin pumping of current. Barnes 
and Maekawa^ show that, in addition to being able to 
generate a spin transfer torque,^ a spin-dependent Berry 
phase can generate what we have called spin pumping. 
Their approach also gives a natural way to understand 
the associated emf-to-field ratio, 2^s/e. 

On the other hand, using irreversible thermodynamics, 
it is clear from the time-reversal properties of the thermo- 
dynamic fluxes that both the spin transfer torque and the 
spin pumping emf correspond to irreversible processes. 
The irreversible nature of these quantities is not clear 
from Berry- phase approaches, where one employs cur- 
rents (fluxes) as primary variables, whereas in consider- 
ing experimental quantities, which correspond to thermal 
averaging, these terms must be considered to be driven 
by thermodynamic forces. 
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